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1.0 INTRODUCTION 


In this report we summarize major research accomplishments that have 

been achieved under support of NASA Grant No. NSG 5341 during the first year 

♦ 

of the grant. The research has concentrated in the following areas t 

i) an examination of observational requirements for predicting con- 
vective storm development and intensity as suggested by recent 
numerical experiments 

ii) interpretation of recent 3D numerical experiments with regard to 
the relationship between O' ershooting tops and surface wind gusts 

iii) the development of software for emulating satellite- inferred cloud 
properties using 3D cloud model-predicted data 

iv) the development of a conceptual/ semi-quantitative jnodel of eastward 
propagating, mesoscale convective complexes forming to the lee of 
the Rocky Mountains. 

2.0 OBSERVATIONAL REQUIREMENTS FOR PREDICTING CONVECTIVE STORM DEVELOPMENT 
AND INTENSITY AS SUGGESTED BY RECENT NUMERICAL EXPERIMENTS 

In the December 1979 Interim Progress Report on this Grant, we concluded, 

based on modeling experience at that time, that the following observations 

are needed for the "precise" prediction of the time and location of convective 

storm development, and the intensity and movement or propagation of convective 

sLorms : 

i) the magnitude of lew*-. level, mesoscale convergence of mass and moist 
static energy 

I i) the vert ical distribution of mcaosen le convergence 
ill) vertical profile;; of the horlr.ent .-.1 wind with special emphasis 

on the lower ont -third of the cloud layer and the temporally varying 


sub -cloud layer 


iv) a definition of mesoscale convergence "anomaly producers". 

Since that time, two sets of numerical experiments have been in pro- 

» 

gress. The first experiment discussed more thoroughly in section 3.0 and 
Appendix 1 involves the simulation of a Florida towering cumuli and its 
dynamic response to seeding. The second experiment described by Cotton, 
Tripoli and Knupp (1980: Appendix 2), involves the simulation of the de- 
velopment and organization of a quasi-steady thunderstorm over South Park, 
Colorado. The mesoscale features affecting the formation of this storm 
have been described by George (1979) and the detailed morphology of the 
storm has been described by Knupp (1980) . The results of these numerical 
experiments support the conclusions derived from earlier experiments and 
add new insight into environmental factors contributing to convective 
storm development, and the intensity and movement or propagation of con- 
vective storms. 

The "dynamic seeding" experiments discussed in section 3.0, suggest 
that the intensity of surface wind gusts and their subsequent impact on 
convective storm propagation are strongly modulated by the vigor and 
thermodynamic consequences of entrainment. The vigor of entrainment, in 
turn, is largely controlled by the strength of vertical shear of the hori- 
zontal wind and by vertical stretching of the updraft as a consequence of 
explosive growth. Also, for a given amount of entrainment, the amount of 
cooling is a function of the dryness of the environmental air at the levels 
of maximum entrainment. 

Some of the environmental factors contributing to explosive growth or 
vertical stretching of an updraft are discussed in section 3.0. Since 
strong shear and dryness of lower, mid-tropospheric air favor vigorous 
entrainment and, fhuslv, strong surface wind gusts, those simulations 
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support the severe weather forecaster's use of minimum G g and shear 
intensity as indices of thunderstorm wind gust intensities. In fact, other 
things being the same (i.e., factors contributing to. explosive growth), one 
would expect that the favored source level of downdraft air would be the 
region where the minimum in 0 e and the maximum in vertical wind shear coin- 
cide. Strong shear thus would favor the initiation and maintenance of 
entrainment at a given level, while a minimum in 0 e would favor strong 
evaporation which initiates and maintains a vigorous downdraft. 

The attempts to simulate the quasi-steady thunderstorm analyzed by 
Knupp (1980) have been quite frustrating. It appears that the scenario 
of prior convective activity over South Park and the horizontal inhomo- 
geneity of the mesoscale convective gust frontal air, all contribute to 
creating the unique organization of the observed storm. A number of 
modelers have speculated that continuously propagating storm systems should 
lose their memory of the cloud initiating mechanism once the main updraft/ 
downdraft couplet becomes well established. In this case, however, the 
observed wind field led to continuous propagation on the N.E. flank; 
whereas it was observed to propagate on its N.W. flank. It appears from 
George's (1979) mesoscale analysis that flow diverted around earlier cells 
could have created an initial N.W. flow at low levels that would have 
initiated propagation on that flank. These results suggest that obser- 
vations of horizontal inhomogenei tics in low’ level flow on the mesoscale 
are needed to uniquely predict a given convective storm, 

3.0 OVERSHOOTING TOPS AND .SURFACE WIND GUSTS OR DOWN BURSTS 

Recently, Cotton, e t a 1 . (1980; Appendix 1) have applied the CSU, 3D 


cloud model to the simulation of the dynamic response of a Fieri Jn cumulus 
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to seeding. ThiB^particular versiqn of the model contains a*, ioe-phaae 
parameterization developed by Stephens (1979). The numerical experiments 
were designed to investigate the cloud dynamic responses to artificially- 
induced "explosive growth". Seeding was accomplished by setting the ice 
crystal concentration at 100*.“ 1 for a period of 10 minutes over the -5° to 
-10°C range and over a horizontal area which constituted the core of the 
rising tower. 

As a consequence of ..-Jing, updraft speeds increased from 15m s“* 
to 33m s -1 and the cloud e>hibited strong overshooting into the stratosphere. 
The cloud penetrated up to 13.9 km MSL and developed strong negative 
buoyancies of up to -6°C at 13 km MSL. Since vigorous tropospheric over- . 
shooting was simulated in this case, it seemed reasonable to expect that 
further analysis of this simulation may shed light on the relationships 
between overshooting and surface wind gusts. One would not expect that 
the intensity of surface wind gusts in the Florida environment to approach 
"downburst" (see Fujita and Byers, 1977), magnitudes, nonetheless the 
mechanisms of communication between tropospheric overshooting and surface 
gusts should be similar. 

In the seeded cloud simulation, a downdraft reaching down to the sur- 
face was predicted, with maximum values of 10.5m s -1 at a height of 1500m 
and 5.9m s'" 1 at 750m, and an associated maximum horizontal divergence at 
0.4 km of 1.24 x 10 2 s _1 . Not surprisingly, this does not quite meet 
Fujita and Byer's (1977) definition of a downburst as "a localised intense 
downdraft witli vertical currents exceeding a downward speed of 12 ft s _1 
(^ 4m s -1 ) at 300 ft O'* 100m) above the surlace." This value corresponds 
to a point divergence of 4 x 10 -2 s" 1 . However, these downdraft velocities 
are significant enough to warrant a closer look at the dynamic links between 
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overshooting tops and downdrafts. 

In the seeded cloud simulation, the overshooting top first manifests 
itself after 2800s of simulated time. At that time, the cloud possesses 
a negative buoyancy of -1°C at 11.6 km. The top continues to rise until 
it reaches a height of 13.9 km at 3000s. It then collapses back, as ex- 
emplified by a maximum downdraft value of 2.3m s -1 at 13.9 km, 100s later. 

A peak downdraft magnitude of 9.9m s -1 in the upper troposphere occurs at 
3400s at an altitude of 10.1 km. This is considerably less than the 41m s 
sinking speeds of the collapsing tops reported by Fujita (1974), however. 

Fujita and Byers (1977) hypothesized the following cell model to ex- 
plain the dynamic linkage between overshooting tops and surface wind gusts 
(dowubursts) . The model involves tops overshooting the anvil then col- 
lapsing into a strong downdraft (located at approximately 10 km) and trail 
of precipitation. Entrainment at the top transports dry air and large 
horizontal momentum downward. Ice crystals are hypothesized to rapidly 
sublime in the subsaturated entrained air, thereby taking up heat from the 
air resulting in a cold, negatively buoyant downdraft. The collapsing top 
and entrained air accelerate the train of precipitation and import fast 
horizontal momentum from the stratosphere. A successive rise and fall of 
the top will create a family of downburst cells that moves away from the 
parent thunderstorm. 

In our simulation, it is not possible to explain the low-level down- 
draft at 1500m, which exceeds 10m s -1 at 30C0s, as the extension of the 
downward motion at cloud top. A very crude trajectory estimate employing 
maximum downdraft velocities for a given height and time, places an upper 
limit of about 5 km for the source level of the sub-cloud downdraft air. 
The bulk of Lhe downdraft air originates from levels substantially below 
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that, as can be seen from cross sections of the flow field. Both numerical 
simulations (Miller, 1978; Tripoli and Cotton, 1980) and observational studies 
(Kropfli and Miller, 1976; Brandes, 1977) have shown that sub-cloud con- 
vective downdrafts originate at midlevels, and not at the tropopause level. 

It thus seems that we can consider the results of this study to be of more 
general validity. 

The fact that in the "no-seed" simulation both sub-cloud downdrafts 
and overshooting tops were much less pronounced than in the "seed" simula- 
tion suggests there is a linkage between overshooting tops and surface 
wind gusts (downbursts) , although different from that envisaged by I’ujita 
and Byers. Based on the dynamic differences between the simulated seeded 
and nonseeded clouds, the following model is proposed. 

Associated with the additional release of latent heat in the seeded 
cloud region is a buoyancy pulse which leads to "explosive" vertical growth. 
One consequence of this increased vertical momentum is the strong over- 
shooting into the stratosphere (see Fig. 1). Another consequence of the 
"explosive" tower growth is the vertical mass flux divergence at and below 
the seeding level, which lends to dynamic entrainment of dry environmental 
air. The strong dynamic entrainment is clearly evident in Fig. 1 as well 
as in the horizontal cross section at 4.9 km shown in Fig. 2. As evidenced 
by the cut off tower, the model appears to overpredict the rate cl dynamic 
entrainment. This is probably a consequence of the design of the seeding 
routine or to the fact that the present model formulation requires that a 
grid volume be wholly saturated or unsaturated whereas the entrained region 
can be expected to be composed of fluctuating saturated and unsatnrated 


blobs of air. 
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As a consequence of the dynamic entrainment of dry environmental air, 
cvaporational cooling, together with the load of liquid water, lead to the 
formation of a downdraft. The downdraft subsequently descends to the sur- 
face being continually reinforced by evaporation of cloud droplets and pre- 
cipitation and by water loading. Thus, the basic mechanism linking over- 
shooting tops and surface wind gusts is dynamic entrainment associated with 
vertical mass flux divergence due to "explosive growth". 

One is led to ask, therefore, what are the factors that naturally 
lead to "explosive growth" and its dynamic consequences? The modeling ex- 
periments briefly described herein, as well as those reported by Tripoli 
and Cotton (1980) suggest that explosive cloud development is strongly 

favored by the presence of mesoscale convergence and by anomalies in sub- 
cloud vertical motion on the cloud-scale. It is interesting to note that 
the case reported by Fujita and Byers involved the partiepitation of a 
sea breeze convergence line as well as an intersecting arc cloud (gust 
front) which provides both the necessary mesoscale and cloud-scale forcing 
mechanisms. 

The numerical experiments reported by Cotton (1972) also suggest that 
the microphysical factors favoring explosive growth would be a suppression 
of a warm cloud precipitation process (collision and coalescence) until 
a tower penetrates deep into the supercooled layer. Tills would allow the 
vertical transport of large qu ntities of supercooled water which is then 
available for subsequent freezing and latent heat release. If the liquid- 
liquid collision and coalescence process then commences, a volatile mixture 
of supercooled raindrops and an exponentially increasing concentration of 
ic.c crystals (as the tower rises and cools in the supercooled region) ran 
lead to the rapid and complete glaciation of the rising tower. The 


V. 
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subsequent release of latent heat would lead to explosive growth In a 
manner analogous to our seeding simulation. 

There are several ramifications of this model of the linkage between 
overshooting tops and surface wind gusts. First of all, from a now-casting 
point of view, the satellite detection of an overshooting top would provide 
little lead time, to the dcsemination of surface wind gust or downburst 
warnings. This is because the formation of the overshooting top would 
occur in concert with the formation of a vigorous downdraft in the lower 
troposphere. Secondly, to forecast downbursts or several wind gusts the 
emphasis should shift from interpreting upper tropospheric environmental 
properties favorable for downdraft formation to lower, mid-tropospheric 
properties. Clearly, strong vertical shear of the horizontal wind in the 
lower, mid- troposphere would favor vigorous entrainment. Furthermore, low 
values of 0 e in the same region would favor vigorous evaporation of cloudy 
air in the entrained regions. 

Clearly .Mils work represents only a small step f orward in our under- 
standing of the factors contributing to explosive growth and the linkages 
between explosive growth in the upper troposphere and surface wind gusts 
or downbursts. It would certainly be desirable to extend this research 
by simulating explosive growth in an actual downhurst-produclng environment. 

4.0 Til Jv DEVELOPMENT OF SOFTWARE FOR EMULATING SATELLITE- INFERRED CLOUD 

PROPERTIES USING 3D CLOUD MODEL- PREDICTED DATA 

Before leaving our group and joining O.E. MATSCO to work with NASA, 

Mr. Mirk Stephens began the development of algorithms lo emulate the analy- 
sis of satellite-derived cloud top isotherms using data predicted with the 
31) cloud model. The algorithm involves n search Cor cloud top and mapping 
of a satellite plane of data (t.c,. 


cloud top temporal are) . Using data 


predicted by the model In its non-Icephase version, he then compulud the 
rate of expansion of isotherm area (— ——) • 

Stephens then compared the predicted — to the predicted hori- 
zontal velocity divergence D. The divergence was computed at the height 

having an average horizontal temperature corresponding to the Isotherm 

d In A 


used in calculating 


dt 


Also, the D was evaluated over the area of 


the isotherm projected onto the plane on which D is to be computed. 

The rate of isotherm expansion was also compared with predicted total 
eddy kinetic energy and the location and amplitude of maximum updraft 
speed as well as precipitation intensity. 

In the test case, Stephen's found generally poor correspondence be- 

J 1 _ » 

tveen — r- — and D. The best correlation occurred at the rime tha: the 

dt 

storm was at its peak kinetic energy early in the cloud's life history 
(30 to 40 rain). After this time, other factors affecting isotherm expan- 
sion became more influential. Seme of these factors were mouel dependent, 
especially for the stage of model development at that time. Factors af- 
fecting isotherm expansion other than velocity divergence Include: 

1) The anvil not only expands outward from a central point « but also 
extends its area by cloud water condensing or moving upward through 
a region. This effect can be expected to be of greater importance 
when the ice-phase is included in the. model predictions. 

2) Since the anvil is composed of cloud water, it is more influenced 
by evaporation than a similar anvil composed of ice particles. 

3) The nhovc-ment toned evaporri *on of Ihe anvil also leads to a 

disorganized anvil structure which contributes to a noisy cr.ti- 
, d In A 

mate of — • 

dt 

The calculation of - - ~ v ~— must also be terminated when the auvij 


dt 


L 

is OUMJO 




I 



reaches the lateral boundary of the model. 

Mr. Dave Hahn (he has recently completed requirements for an M.S. 
degree) ia now continuing the development of algorithms emulating satellite 
views of convective storms. This work will make use of model-predicted 
data in which the ice phase is included in the model (sec sections 2.0 and 
3.0). 



5.0 DEVELOPMENT OF A CONCEPTUAL / S EMT - QUANT I TATI VE MODEL OF EASTWARD PRO- 
PAGATING, MESOSCALE CONVECTIVE COMPLEXES FORMING TO THE LEE OF THE 

ROCKY MOUNTAINS 

The development of a conceptual/scmi-quantitative model of eastward 
propagating, mesoscale convective complexes has been concentrated on the 
analysis of convective complexes forming during the period 3 August to 
10 August, 1977. A two-part paper is now being prepared for submission to 
the Monthly Weathe r Review entitled: A Long-Lived Mesoscale Convective 
Complex, Part I - The Mountain Generated Component, and Part II - 
Morphology of the Mature Complex. 

Part T is basically a synthesis of George's (1979) thesis. The paper 
follows a scenario of convective envolution from the early morning micro- 
meteorological scale, to the raid-day ridge/valley circulation raid Mountain 
v umull /cumulonimbus , to mid-afternoon plains cumulonimbus and squall lines, 
to the formation of a mesoscale convective complex (MCC). 

In Part II, the pre-MCC synoptic field is re- analysed and the morphology 
of the MCC is described including: 

i) the interaction of the MCC with tins reunins of the previous day's 
MCC 

ii) the intrusion of dry, high momentum air into the back ol' the storm 
at. 500 mb 
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ill) a broad region of high momentum outflow from the storm just below 
the tropopause 

iv) a significant raeso-anticyclone at 200 mb which was not present prior 
to the formation of the MCC 

v) lack of dire.ct interaction with the polar front jet well to the 
north 

vi) divergence profiles similar to tropical clusters 

vii) virtually no evidence of poleward meridional heat transport asso- 
ciated with the MCC as one would expect in a baroclinic system. 

It is concluded that the MCC's studied are basically tropical in 
nature and their dynamics are dominated by buoyant accelerations. Asso- 
ciated weak, baroclinic features, such as a surface stationary front, serve 
to trigger and direct the release of convective instability. It is sug- 
gested from this analysis that a certain amount of baroclinicity may be 
compatible with an MCC, but as the baroclincity increases the convection 
either tends to organize itself into the linear structure of a squall 
line or the system as a whole undergoes a transformation into a rapidly 
occluding cyclonic wave. 

Dr. Pete Wetzel is a major contributor to Part II of these papers, 
therefore our goal is to complete these manuscripts before Dr. Wetzel joins 
the NASA/GLAS Severe Storms Group in Sept. 1980. 
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1. Introduction 

"Dynamic seeding" of convective clouds as it is being performed in the NOAA Florida 
Area Cumulus Experiment (FACE) rests on the assumption that seeding of actively growing 
tc.-ers will lead to a conversion of supercooled liquid water to ice and subsequent 
"explosive growth" due to the latent heat release. The accelerated growth of rising 
towers is then hypothesised to lead to strengthened low level inflow, subcloud layer con- 
vergence and possibly even merger of neighboring clouds. It was pointed out by Cotton 
and Tripoli (1979) and Simpson (1980) that downdrafts might play a key role in communi- 
cating enhanced buoyancy at the -10°C level down to the subcloud layer. 

The present study will attempt to provide answers to some of the questions raised 
with respect to the dynamic response of a cumulus to seeding. A three-dimensional model 
which includes an ice-phase ’parameterization is used to determine differences between a 
seeded and unseeded cloud. In this paper, a short summary of the model is given, the case 
study is described and the results of the numerical experiments are discussed. 


2. Model Description 

The three-dimensional cloud model employed in this study is a revised version of the 
model i .ported by Cotton and Tripoli (1978) and Tripoli and Cotton (1980). Major changes 
include the use of the quasi-conservative thermodynamic variable ice-water potential 
temperature 0^ (Tripoli and Cotton, 1979) and the addition of an ice-phase paraneterizatie: 
(Stephens, 1979). Some other major features of the model are: 

1) it is fully 3D 

2) it is non-hydrostatic, time-split compressible 

3) it utilizes sigma-z vertical coordinate 

4) it presently has "Smagor insky-type" eddy viscosity, first-order turbulence closure 

5) its lower boundary condition on turbulent fluxes is derived from a similarity 

theory based surface layer parameterization 

6) it contains a simple rain parameterization using autoconversion and accretion 

assuming rain is distributed according to: 


•I’ (10 = “ I •■XI 1 (-K/\) 
m 


with R constant 
m 


. , * 

- ■' • V 

■..•■■■••a 




it contains an ice-phase parameterization developed by Stephens (1979) 
the lower vertical boundary condition is given by a surface layer turbulence 
parameterization and zero mean advective fluxes. The top boundary is a 3 kn 
deep Rayleigh friction layer capped by a rigid lid 
the lateral boundary conditions are a variation on the Klemp-Wilhelmson (1978) 
radiative lateral boundary condition 
it uses second-order finite differencing in space 

its time differencing scheme is Matsuno (1966) on small time-steps and leap-frog 
with Asselin filter on long time-steps' 
its grid-spacing is 750 meters in the horizontal and vertical over a 36.75 x 
36.75 km (horizontal) x 15 km (vertical) domain 
it is initialized with a base-state sounding, large-scale convergence and a 
dynamic focusing algorithm described by Tripoli and Cotton (1980). 


3. Case Study Description and Design of Experiment 

The numerical experiments are performed in an environment described by measurements 
taken on the 25th of August 1975 as part of NOAA/Florida Area Cumulus Experiment (FACE) . 

Two convective systems formed inside the surface mesonetwork on that day. The second ens 
of these was an isolated system and is subject of a case study by Cunning et al . (1979). 

It formed at about 1S30 GMT and three of its growing towers were seeded with silver iodide 
(Agl) by aircraft at 1834 GMT, 1838 GMT and 1S40 GMT. Each tower was seeded when it 
reached the -10°C level. Following seeding, a fourth (unseeded) tower grew explosively to 
a higher level chan either one of the three previous cells. Throughout its life-cycle, 
the system was continuously scanned by triple Doppler radar. 

The present study will not attempt to reproduce the features of any particular cloud, 
but rather investigate the effects of seeding on a cloud of comparable dimensions in a 
similar environment. Accordingly, the base state sounding used to initialize the model 
is a modified 1745 GMT special sounding taken at the center of the observation network 
matched with the Miami 1200 GMT sounding above 300 mb. The wind profile is composited 
from low level PIBAL observations at the central observation site at 1745 GMT and Miami 
1200 GMT rawinsonde data above that. To prevent the cloud from being advected out of the 
domain, the low-level mean was removed. The sounding is conditionally unstable and 
relatively dry above 700 mb. The vertical wind shear which is most pronounced in the 
u-component can be seen from Figure 1. The imposed mean vertical motion field is based on 
numerical model rr Its by Tielke (personal communication) and is matched to an observed 
value of surface convergence of 2 x 10 - ^ s“^. Numerical values are 0.09 m s“l at 750 m, 

0.13 m s -1 at 1500 m, 0.10 m s“l at 2250 m, 0.05 m s~l at 3000 m and 0.02 m s“l at 3750 m. 
Twenty-eight percent of this convergence was focused on a scale of 2.25 km, resulting in a 
peak vertical velocity of 3 m s - * in the center of the domain. Experiment I was initialize, 
with this procedure. Experiment II. is identical to Experiment I except that seeding is 
simulated after 41 min 40 sec, at which time the cloud has readied the -10°C level. Based 
on the number of flares released into the cloud system and their nucleation efficiency 
(Sax et al ., 1979) the total number of artificial nuclei was estimated. In consistency 
with this estimate, the ice-nuclei concentration was set equal to 100 £“1 over a peribd of 
10 min for 2 levels in the vertical (~5°C to -10°C range) over a horizontal area which 
constituted the core of the rising tower (3.75 km < x < 5.25 km, and -1.5 km £ y <^-.75 km) . 

ORIGINAL PAGE IS 

4. Results OF POOR QUALITY 

Experiment I: 


Due to the pure convergence initialization without any temperature or moisture rer- 
turbat ion the initial growth ol the cloud is fairly slow. At 30 min, the cloud reached 
a cloud top height of 4 km, after which time the cloud top rises hv about 5 n s" ^ . The 
vertical wind shear leads to a Lilted updraft and strong entrainment of dry air on the 



Fig. 1: Cloud development at a) 45 minutes and b) 50 minutes. The thick solid line 

narks the cloud outline, thin solid lines are contours of total condensate (every 2 g kg - ^). 
The presence of rainwater and graupel is denoted by a dot or a triangle, respectively. 


dovnshear side of the cloud (Fig. 1). After one hour, the rising tower is completely 
cut off from its low-level support and rises in a bubble-like fashion up to the tropo- 
pause. The "t ime-window" effect for seeding potential is evident from the fact that 
virtually all the supercooled liquid water at the -10°C level is frozen within a time 
period of less than 10 min. In comparison with the observations, it should be noted that 
due to the limited resolution of the model, only one large rising tower is reproduced 
rather than four individual cells. 

Experiment II: 




"Seeding" the rising tower not only resulted in a complete glaciation of the cloud 
throughout the seeded volume mainly by conversion of liquid water to graupel, hut produced 
marked dynamic changes as well. Updraft speeds increased from 15 m s“l to 33 m s”l and 
the updraft widened and deepened. The correlation between vertical explosive growth and 
widening of the cloud was documented by Simpson, ejt_nJ_. (1963) and Simpson, Brier and 
Simpson (1967). The pressure minimum at middle levels associated with the updraft was 
much more pronounced following seeding while at the same time pressure increased slightly 
below cloud base. The downdraft on the down she nr side of the cloud was increased in 
intensity. These differences in the pressure field and draft structure rem..ined very 
pronounced over 10 min following seeding. At that time, the tower which consisted mainly 
of ice crystals had a top about 3 km higher than in Experiment I. At one hour simulated 
time ( LS min after seeding) the tower again was compl etely separated from the low levels 
and the intensity of the circulation had fallen hack to the levels of Experiment 1. 

Total prec ip Ltat ion was increased by almost 100E over Experiment T. 
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Conclusions 


It should be noted that the renults presented in this paper are of a preliminary 
nature and subject to revision after a more detailed analysis of the numerical experi- 
ments. At this stage, the main effect of seeding is found to be accelerated tower gro 
and an increase in precipitation, while the link to the subcloud layer remains unclear 
The intensification of the downdrafts failed to initiate or sustain convection in this 
case study. The reasons for this are currently being reviewed and will be presented a 
the conference. 
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Observational Analysis of an Intense, Quasi- 
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1. Introduction 

On 19 July, 1977, an Intense, quasi-steady 
thunderstorn was observed by throe Doppler ra- 
dars, ravinsendes, and tlm NCAR Portable Auto- 
mated Mcsonct (PAM) over Soutn Para. Colorado, 
during the Colorado State University. 1977 
South Park Area Cumulus Experiment (SPACE-77). 

In this paper, we summarize the salient 
mcsoscalc featurea leading to the foimacion of 
the quasi-steady thunderstorm and the observed 
structure of the storm. The results of numeri- 
cal experiments with the three-dimensional cloud 
model reported by Cotton and Tripoli (1978), 
Tripoli and Cotton (1980) and Cotton, et nl. 
(1980) attempting to simulate the formation 
and structure of the observed storm are 
discussed. 

2. Observations 

2.1 Mesoscale environment and general echo 
characteristics 

The 19 July regional synoptic environment 
was characterized by relatively weak-shear 
5-10 m s"* southerly flow at mid to high levels. 
Throughout the day, the NE transport of low- 
level moisture over Colorado resulted in early 
deep convection over NW Colorado and nidafter- 
noon deep convection over South Park. An analy- 
sis of the evolution of the 19 July South Park 
mesoscale features by George (1979) revealed a 
relatively complex, time-evolving pattern. 

During the initial stages, a norrh-couth line 
of echoes was Induced by a corresponding north- 
south line of low-level convergence, with 5-7 
m s' 1 westerly winds to the west, and 5-7 m s'l 
easterly winds to the ease. New convertive 
cells 6-8 km In diameter typically formed 
periodically (10-20 min) on the line’s southern 
(upsiiear) end and moved northward at ft-R m s"l 
through the line while slowly weakening. 

Hie low-level flow and echo patterns at 
1742 MDT (approximately 90 nln attcr the line's 
formation) are portrayed In rig. I. At this 
time, alr.'iow behind a southward- ad vanning 
neso-cold front was characterized hv relatively 
strong northerly surface. winds (vnpeci.it lv to 
the west oi the line). The vertical structure 
of this foul— frontal air mass exhibited sic.nl- 
f leant lv Tenter low-level shear at the hoi i- 
gunlil wmJ i .|il*J s* * ) and relatively I area 
1 -w- level “ixu.e ratios ( Mg kg"*). I he totJl 

huovanev ,, |>jreei lilted iron the start ace 
.as similar on each side o| the : r 'nt . iliiiough 
reater ! w- level negittve nuevanev existed In 


the cooler air. The post-frontal northerly flaw 
remained over the western half of South Parkier 
the er uing 90 min. The increased low-ievel 
shear of the relatively shallow (£1 km) airmass 
apparently altered subsequent storm structure 
and motion. One preexisting cell within the 
echo line's Interior spawned s secondary ceil 
which split and exhibited a diverging trajec- 
tory. Another storm, discussed in greater de- 
tail in the following section, rapidly intensi- 
fied upon encountering the moist, northerly flow. 

2.2 Observed characteristics of a quasi-sleadv 

storm (till) 

The most significant chance in post- 
frontal storm behavior was the rapid intensifi- 
cation and organization of a relatively weak, 
multicellular cluster located on the southern 
end of the echo line. Prior to the arrival of 
the neso-cold front, this cell group exhibited 
a transient behavior in echo characteristics. 
Contrastingly, after the front's passage, the 
cell cluster rapid.’y intensified and consoli- 
dated within 30 rain to become a heavy precipi- 
tating, quasi-steady storm which traveled \i0 %> 
to the left of the mean environmental cloud- 
level winds for the next 60 rain. 





Fig. 1 Surface t !**w streamlines, nenonet 
parameters, and ».V i'P- 1 LPl echo contours iJ5 
.iliZ and Art .132) at 17AJ ill'T. Tr.'l 1i clones •'! 
echo centroids arc depleted. The 1710 vertical 
sounding lee it Ion is denoted l>v star. (From 
George, Id*’)). 
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The mil lent feature* of the lieu pattern* 
within and adjacent to Cll during Ita Hteadv 
period are portrayed In Fie. 2. Storm notion 
was directed toward* the KU in response to con- 
tinuous regeneration of precipitation within 
the propagating updrart Ui, located above tho 
gust front in the NW storm quadrant. The 
steadiness of storm circulation patterns is at- . 
trlbuteu to the constructive Interaction of the 
KU-euvinq gust front with opposing low- level 
flow. The observed behavior is consistent with 
previous numerical experiments (eg., Klemp and 
Wllhelrason. 1978; Thorpe and Miller. 197S) which 
have elucidated the inporcanco of wind shear 
profiles on storm oat ion and organization. Air- 
flow at mldlevela was characterized by acceler- 
ated (low around the updrafts and a weak flow 
vaka region extending downstream. Anticyclonlc 
vorticlty puaking at mldlevela (1.2 x 10'2 * -i ) 
was associated with updraft Ul, while cyclonic 
vorticlty was associated with flow around the 
updraft U2. Tilting of vortex tubes was a 
primary contributor to these patterns. High 
turbulence (Doppler radar Inferred) was es- 
pecially pronounced and persistent at nidluvrls 
along the southern and western storm quadrants. 
Estimated turbulent kinetic energy dissipation 
rates exceeded 0.1 m2 s“3 within a highly 
sheared region separating an updraft and down- 
draft. Less intense turbulence was associated 
with mldlevel relative storm inflow along the 
southern storm quadrant. 

Doppler-derived draft patterns (Fig. 2) 
exhibited more complexity than corrr»pwndlng 
reflectivity patterns. Updraft Ul had roots 
ahead of the gust front in the NW (downxhear) 
and exhibited a reversal in tilt from lower to 
upper levcle. Continuous northwestward pro- 
pagation of Ul apparently goverened ntorii move- 
ment. Analyzed maximum updraft speeds attained 
peak magnitudes of 20-30 a s~l at upper mid- 
levels (708 knACL). A secondary, less intesnse 
and extensive updraft, with peak speeds of 10- 
20 m s"f , was analyzed in the southern storm 
quadrant. Its forcing mechanism could not be 
resolved from the Doppler data or surface meso- 
net data. Diverging of low-level flow around 
the west side of the gust front and subsequent 
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Fig. 2 Concept u »l si'de 1 of Cll's :l»w pat- 
terns during its intense, organized don stage. 


convergence along the southern flank, ainilar to 
what has been modeled by Tripoli and Cotton 
(1980) may have plaved a role. The analyzed 
downdraft circulations consisted of three cells 
or sourca regions which coalesced near the sur- 
f«c« to product a coherent cold air mass and 
associated gust front. Inferred downdraft 
Initiating and maintenance mecnanism* consisted 
of precipitation drag (Dt) and evnporational 
cooling (02 and Dl). Flow around U2 waa instru- 
mental in producing downdraft D2. 

One of the most striking characteristics 
o( Cll waa ita copious production of precipita- 
tion. Of tha 12 primary storm cells observed 
during a 44-hour period within South Park. Cll 
contributed •OOZ of the total precipitation. 

Such a large fraction waa attributed to Cll'e 
quasi-steady organized circulation and relative- 
ly long lifetime. The peak measured rainfall 
rate (1833 MDT) was 190 mm hr -1 , corresponding 
to a maximum low-level reflectivity of less 
than 30 dBZ. Even higher rainfall rates may 
have existed after 1900 MDT since lew-level re- 
flectivity factors exhibited an Increasing 
trend until 1937. Analysis of Cll's radar echo 
characteristics revealed that areas of 40 dSZ 
and 43 dBZ echoes at 1 km ACL, as well as vol- 
umes enclosed by 40 dBZ and 45 dBZ echo surfaces 
attained maximum values at 1930. In contrast, 
updraft (and mass flux) magnitudes decreased 
slowly after attaining peak strength at 1900 MDT. 
Such a relationship implies that Cll's precipi- 
tation efficiency increased with time. 

.Despite Cll's dynamic vigor and organisa- 
tion, large hall was inferred absent because 
of the lack of significantly hign reflectivity 
(53 dBZ maximum). However, small hall znd/or 
graupel may have reached the surface. The ab- 
sence of large hall, together with the presence 
of high rainfall rates. Implies that precipita- 
tion processes within Cll were extremely effi- 
cient. Analyzed airflow and reflectivity pat- 
terns suggest that such a high precipitation 
efficiency resulted from an organized recircu- 
lation of precipitation from the updraft at aid- 
levels to that at low levels. Recirculating 
tranjcctorlca were made possible bv the large 
angle (n.140 0 ) between the mldlevel and low- level 
wind vectors. Therefore, ait flow and turbulent 
notions adjacent to the midlevel updral't would 
have transported precipitation from the updraft 
northward over the low-level inflow. The high 
precipitation efficiency also implies that the 
liquid water depletion concept advanced bv 
Foutu (1979) was occurring naturally with Cll, 
where numerous particles reentered the down shear 
updraft (Ul) in both liquid and solid fern (tfcii 
la substantiated by visual observations ot a ne- 
bulous inflow sector). Thus, recirculation of 
piecipitatlen elements explains (a) high rainfall 
rates, (b) lack ot large hail, and (c) absence 
of a radar echo vault. 

In summary, the change in luw-|evel air- 
flow affected storm-scale processes in two wavs: 
(1) It altered torn dynamical processes nv 
producing a quasi -steady iipdr.il t which governed 
storm movement; i.') the low to nldlevel level 
directional wind shear pr«mot ed in orcinlzeJ re- 
circ'ilatlen oi precipit.it ion elements into tne 
l.-w-level updralt, thus nod iivuu precipitation 
growth processes. 




3.0 Numeric ;il Experimental 

In order to yciiii j greater understanding 
of the dynamics responsible for lira product ton 
and maintenance of Ell, a three— dimensional 
numerical simulation of cumulus development in 
the observed 17 July, 1977 environment was per- 
formed. The Improved CSIi multidimensional cloud 
model summarized oy Cotton et at, i i 9 SO 1 was 
used. This model utilizes an ice phase para- 
mcterlzat ion developed by Stcpnens il9'9) in 
which mixing ratios of lcc crystals and larger 
qriupei particles are explicit!*-* predicted. 

Because the ultimate objective of this 
study is to use cue modal to gain further in- 
sight Into the dynamical structure of C 1 1 . the 
basic state temperature, moisture end wind 
fields were specified initially from a late 
afternoon post-frontal sounding. A domain 
of J5 km square by 17 km high was selected with 
a uniform 750 n grid spacing. A ronvective cir- 
culation was initiated by the Introduction of 
a low-level focused convergence field (see 
Tripoli and Colton. 1980) producing a strong 
Induced uplifting within an A km radius of lira 
domain center. Surface moisture and potential 
temperatures were iniciallv perturbed upward 
to the level of free convection in ihe region 
of maximum convergence so that downward accelera- 
tion produced hy negative buoyancy would not 
destroy the perturoation before the cloud became 
self-sustaining. As a result, an intense cumu- 
lus circulation developed reaching 14.5 kia MSL. 
The updraft, with a peak mignituHe of 14 D s*l, 
evolved from the. initial upward motion and mi- 
grated to the southern quadrant oi the storm. 
Tills migration can be explained hy the findings 
o£ Tripoli and Cotton (I960; hereafter referred 
to as T-C), which showed that in the presence 
of strong Initial low-level convergence, the 
perturbation opar ift may become long-lived. As- 
sociated with viperous well-rooted, initial dis- 
turbances was a surface pressure low. In some 
cases, T-C found that the pressure anomaly was 
Intense enough to divert downdraft air beneath 
the updraft in a direction opposite to the 
downward transport of horlaont.il momentum. Thus 
the downdraft reinforced the initial updraft 
which In this case was located In the southern 
quadrant ot Ihe storm. 

Because the updraft hi C!1 was, in fact, 
observed in the northern quadrant of the storm, 
it was decided to try a different approach to 
storm initialization. Storm i'll was observed 
to first develop in the prefrontal environment 
he for v taking on the ch.ir.ieter tat les of a steady 
state c lrcul.ilion subsequent to front. il Passage. 
Tills earlier development took nl ice along the 
soutiierr flank ol a convective line tarred in 
lonluncll >n with a preexisting nesisrale con- 
vergence zone. In addition, the Jaw-1,- •'•el wind 
aiear wis apposite ta tint encountered l iter on. 
It was decided, there to re. to simulate convec- 
l lan r 1 r -it in tnc preireniai environment an ! 
tacit arth I. tally Induce a tr -ul il passage ;r. 
the presence at mature c.iin'.-it tVe .lilj. 


was again taken, tnlt hilly, a uniform surface 
mesoscale convergence of 0.77 x 10"’’ s“* (simi- 
lar Lo observed) was Imposed, A region of -am 
vigorous vertical notion was imposed In the 
center of the southern 2 ' J of th« domain Vlthir 
a diameter of 8 km. The tecusing function was 
slightly different than that of T-C In order to 
contain all compensating subsidence within a 
radius of 15 km. Instead a: over tne entire et- 
na In . Tills was necessary pecausc the perturtt- 
don Is not centered in the domain. 1: was 
anticipated that with this Initialization mo:-... 
convection will develop in conjunction with :..i 
focusing, drift northward and develop new 
cells tp the south as observed. The result! 
show that Indeed convection developed, drift*! 
northward and some new convection was produce! 
to the south. The initial forced cell reach*! 
14.5 km MSL which can be compared to 15 km Mil 
observed in the prefrontal environment. Peak 
updrafts of 24 m s**l were produced. No pre- 
frontal estimates of vertical velocity were mad* 
from observations, however values of 25 m *“* 
were observed in Cll after frontal passage. Pre- 
cipitation of over 1.5 cm over a 3.5 km- area 
occurred primarily in the form of zrsupcl. This 
is compared to ohserveo values of 1.0 cm, al- 
though sampling stations rav not have been lo- 
cated in the maximum zone. Fig. 3 shows the 
cloud wntor field at 2250 n above the ground 
level at'cer 50 sin simulation time. A loosely 
organized north-south pattern of convective 
cells cun bo seen. 

Numerical experiments in Progress at the 
time of this writing include an altered initiali- 
zation algorithm designed to tighten up the 
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urn. inln.it loti of the north-smith lino oi col In. 
In addition, the passu no of the no --cold 
front through the lino or nature convective 
cells will bo induced. The toitults of these! 
nunoric.il experiments will bo described in 
tlio or.il presentations. 
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